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Abstract 

Rotational spectra of diatomic molecules measured in the high-precision 
experiments are analyzed. Such a spectrum is usually fitted by an 8th 
order polynomial in spin. In fact, from the theoretical point of view, the 
rotational spectrum is expected to be a smooth function of spin. Thus, 
fitting the rotational spectrum in terms of a low order polynomial in spin 
seems to be justified. However, this is not what we found in the present 
work. We show that the measured spectrum is not necessarily a smooth 
function of the spin, which is quite surprising. Whether or not it indicates 
the presence of some hidden physics is an open question. 
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I. INTRODUCTION 



Characteristic quantities of a diatomic molecule are the energy associated with the 
motion of electrons and the one associated with the motion of two nuclei (ionic cores). 
The latter motion is essentially that of a dumbbell shaped system (vibration against one 
another and rotation as a whole). The order of magnitude of the electronic vs. vibrational 
excitation energy {Eeie/ Emb) is typically 100 and that of the vibrational vs. rotational 
one {E^ib/Erot) another 100. It means that {Erot) corresponds to transitions in the 
near (far) infrared. Thus, the electronic motion is very fast compared with that of nuclei, 
so that its wavefunction may be constructed from the instantaneous constellation of two 
nuclei assuming as if nuclei were frozen and thus depends parametrically on the positions 
of two nuclei (Born-Oppenheimer approximation). This is the standard way of decoupling 
the electronic and nuclear degrees of freedom from one another. The treatment of the 
nuclear part then becomes very simple. 

The nuclear radial Schrodinger equation (with k — 1) takes the form 



where r is the distance between two nuclei and M the reduced mass. The potential U (r) 
between two nuclei is a complicated object that involves an integration over the electronic 
wavefunction, which we do not go into here. In practice, it is phenomenologically replaced 
by a Morse potential. The quantum number I is the rotational spin which characterizes 
the angular wavefunction Yjm{(^,(P) and V the so-called vibrational quantum number 
which characterizes the radial wavefunction. We note that, because of the mathematical 
property of this equation, the eigenvalue Ev{I) has to be an analytic function of the 
parameter (spin) I. Since the effective potential 



has a prominent minimum at r = Tq, where Tq is of the order of the molecular size, one 
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V{r) = U{r) + 
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usually expands it around r = ro ( 



dVjr) 



0) and expresses V{r) in a power series of 



dr 



J r=ro 
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r — tq. Consequently, the radial motion can be treated as a vibration having a perturbing 
unharmonicity. The dependence of Ev{I) on / becomes then a power series of /(/ + 1) 
and one usually stops the series at the 8th order in /: 

Ev{I) =Tv + BvI{I + 1) - Dv[I{I + 1)]' + Hv[I{I + 1)]' + Ly[/(/ + 1)]^ (3) 

In recent years, very accurate data of rotational spectra in diatomic molecular bands 
became available thanks to modern experimental techniques [1-8]. What is measured in 
such an experiment is the so-called R- and P-process |^ 

R{I) = Eh{I + 1) - El{I), P{I) = Eh{I - 1) - El{I). (4) 

These quantities are the 'inter-band' AI = 1 (dipole) transition energies between a higher 
{V = H) and a lower {V = L) vibrational band and are usually measured in the unit 
of wave number (cm~^). The coefficients in the formula (|^) have been fitted to such 
experimental data [1-8]. However, one can take a different approach ||1CI|| . 

By inverting the relation (^), one obtains the 'intra-band' A/ = 2 (quadrupole) tran- 
sition energy AEv{I) = Ev{I) - Ev{I - 2) as 

AEh{I) = R{I - 1) - P(/ - 1), AEl{I) = R{I - 2) - P{I) (5) 

for the higher and lower band, respectively. It should be remarked that, from the set of 
data of R- and P-process, one obtains four separate sets of data AEy{I) ioi V = H and 
V = L bands with even as well as odd spin sequences. The intra-band transition energy 
should be a smooth function of spin / according to the theoretical consideration made 
in the beginning. In the present work, we want to examine this statement. We will in 
fact show that the measured intra-band transition energies are not necessarily smooth 
functions of spin. Let us first devise a tool which suits the analysis of data. 

II. THE METHOD OF ANALYSIS 

To study the behavior of AI = 2 transition energy AE{I) = E{I) — E{I — 2) for 
each band {V = H or L) and each spin sequence (/ =even or odd), some kind of data 
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manipulation is necessary. In fact, if one plots the quantity AE{I) directly, it is hardly 
possible to see any fine structure since AE{I) is a globally increasing function of / and 
extends over a wide range of values. Thus, the basic idea is to look at the deviation of 
the transition energy AE{I) from its smoothly increasing part. 

To this purpose, we subtract a polynomial of order iV in / from AE{I) and define what 
we call the A^th order one-point formula |ll| 



N 

aS^e(/) = AEii) - g^(/), g^(/) = J2 Qmi"" (6) 

where the coefficients Qm are determined by minimizing the quantity 

x{qo,---,qN)^[A^N^E{I)f (7) 
with respect to g.^ = 0). This leads to a set of + 1 equations (m = 0, 1, ■ ■ ■ , A^) 

N AI=2 AI=2 

J2Smnqn = Tm, Smn = E ^™ = E ^'"A^W- (8) 

n=0 / / 

We note that the smooth part Qn{I) which we subtract from AE{I) is nothing other 
than a polynomial of order determined by the x-square fit to AE{I). However, in 
practice, this formula cannot be used in the present form particularly when the order of 
the polynomial A^ is larger than 3 since the equation (j^) is highly ill-conditioned. Thus, 
in what follows, we want to transform it into another form. 

First, let us note that the replacement I ^ I — Iq does not change the fitting procedure 
since the shape of the polynomial Qn{I) as a function of / is unchanged. Thus, the origin 
of the spin values can be shifted freely. Secondly, the spin variable can be scaled too 
(/ — > al) since the order of the polynomial remains the same. These properties can be 
used to rewrite the polynomial in a different form. 

Shifting and scaling the spin values can be achieved most generally by a linear mapping 
I = ax + b. The increment of x is thus Ax = ^ (A/ = 2). We will choose a = ^'"<^^-^^i" 
and b = so that the range of x becomes [-!,+!], where x = —I (+1) corresponds 

to / = Imin {Imax}- The polynomial in question may thus be written in the form 

N 

Qn{I) = E PmPmix). (9) 
m=0 



Here, we use the Legendre polynomial Pm{x) instead of x"^. The reason will be explained 
below. The resulting set of equations is similar to (|]) but g„ is replaced by Pn and J" (J'") 
by Pn{x) {Pm{x)). This representation has an advantage that there holds the relation 

+1 

Smn = X! Pm{x)Pn{x) = if m + n = odd. (10) 

x=-l 

It means that the whole set of equations splits into two independent sets of equations of 
smaller dimensions, one with m, n =even and the other with m, n =odd: 

N +1 +1 

SmnPn = Tm, Smn = Pm{x)Pn{x), = ^ Pm{x)l^E{I). (11) 

n=even or odd x=—l x=—l 

This set of equations determines the coefficients Pm and accordingly the polynomial (P) 
which represents the smooth part of AE{I). 

It remains to be remarked that the original set of equations (§) is highly ill-conditioned. 
It is indeed so ill-conditioned that even the double precision algorithm is not free from 
the numerical instability caused by large losses of accuracy if is greater than 3. This 
problem can be avoided by shifting and scaling the spin values as presented above. In 
fact, the situation improves slightly if one uses the power series x"^ thanks to a property 
analogous to (pil|). Nevertheless, this does not fully resolve the numerical instability. The 
reason lies basically in the fact that the zeros of are multiple and are all concentrated 
at X = 0. In contrast, all zeros of Pm{x) are simple and never coincide with one another 
for different m's. The use of the Legendre polynomial Pm{x) instead of the power series 
x"^ is thus essential for the numerical reliability (both accuracy and stability). A similar 
caution may be necessary when fitting the formula (^ to experimental data. 

In the next section, we will present A^^^ E{I) for a number of experimental data using 
a 7th order polynomial. In fact, = 7 is quite appropriate for the analysis of molecular 
rotational spectra since the rotational spectrum for a given vibrational band is fitted 
(globally) by a formula containing up to the 8th order in spin / (see eq.(0)) so that the 
transition energy AE{I) is a 7th order polynomial. 

We will also define the 'filtered' one-point formula by setting the quantity A^^^ E{I) to 
zero if its absolute value is smaller than or equal to the corresponding error bar. This 
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formula is quite useful in practice. By construction, it shows whether the deviation of the 
transition energy AE{I) from its smooth part is physically significant or not. 

III. ANALYSIS OF EXPERIMENTAL DATA 

As mentioned before, we obtain four sets of transition energies AE{I) (even and odd 
spin sequence for higher and lower band) from the measurement of R{I) and -P(/). Thus, 
four independent sets of data A^^^ E{I) can be created, which we present as diagrams. 

As an example, we plot in Fig. la the 7th order one-point formula A^p E{I) obtained 
from the measured R{I) and P{I) of the 2-2 band of the A^S^ — system of the 

molecule YN [Q. On the other hand. Fig. lb shows the corresponding filtered one-point 
formula. It is 'filtered' in the sense that we take 

A^^^E{I) = if \A^^^E{I)\ < Error Bar (12) 

where Error Bar is twice the error bar of the measurement of R{I) and P{I), see eq.(^. 
Eq. ( p!2D makes easily visible that some deviations of the transition energy from the smooth 
part that occur in Fig. la are physically insignificant as they are zero within the Error Bar. 
Thus, if the rotational spectrum is indeed a smooth function of spin J, which is what 
we want to confirm, the deviation should be zero everywhere when filtered. However, 
Fig. lb clearly shows that there remain some deviations which are not filtered away to zero 
{V = H band). These deviations are beyond the Error Bar and may thus be considered 
as significant. This is rather surprising but is not a single case. We found many such 
cases, most of them being more complicated than this example. It should be emphasized 
that the filtered formula clearly shows at which spins the anomaly occurs. 

We present another example. Figs. 2a and 2b show respectively the one-point formula 
and the corresponding filtered formula applied to the 1-1 band of the C^S+ — X^S"*" 
system of the molecule YD 0. One sees a behavior similar to the above example, though 
it looks slightly complex. As a matter of fact, this kind of results is seen also in the 0-0 
and 1-0 bands of the — X^S+ system of the CrD molecule p. Actually, there are 
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many more complex cases in the data we have analyzed, some of which will be shown 
later. Thus, it seems that there exists anomaly (irregularity) in the measured spectrum 
and the question is where it comes from. 

We have examined the relation between a spectrum and the resulting irregularity which 
appears in the transition energy. Suppose that there is a band crossing in a spectrum. It 
will produce a kink at the crossing point and thus an irregularity in the transition energy. 
While it is unlikely that a band crossing may occur in the molecular cases, it indeed exists 
in the nuclear cases [ill], which leads to the so-called backbending phenomenon. Based 



on this observation, we have studied a model in which a smooth rotational spectrum 
E{I) = AI{I + 1) has a sequence of kinks around a spin / = Iq. Since the smooth 
part of the transition energy AE{I) = E{I) — E{I — 2) = 2A{2I — 1) is a first order 
polynomial in J, it will be sufficient to use the 1st order one-point formula A^i^ E{I) in 
such a demonstrative example. Some of them are presented in Fig. 3. 

The three-kink model in Fig. 3 is of particular interest to the YN molecule Q as it 
reproduces the V = H even spin data shown in Fig. lb while the V = H odd spin data 
corresponds to a special four-kink model. The three-kink spectrum in Fig. 3 was produced 
by shifting 'down' the energy from a smooth curve E{I) = AI{I + 1) at a spin I = Iq while 
the pattern with exactly opposite phase will be obtained by shifting 'up' the energy. In 
general, a spectrum with a given number of kinks can produce several different patterns 
depending on the way how the kinks are created. 

Let us next examine the YD molecule presented in Fig. 2b. The V = H even 
spin data can also be understood in terms of a three-kink model, in which three-fold 
kinks occur at two different places (at / = 18 and 32 in a spectrum of the form E{I) = 
AI{I + 1) — B[I{I + 1)]^). Fig.4 compares the data and such a 'theory'. The experiment 
is reproduced quite well. On the other hand, the V = H odd spin data seems to suggest 
the presence of different types of kinks that occur at three different places. In fact, the 
first one (at the lowest spin region) shows obviously a two-kink pattern while the last one 
(at the highest spin region) a four-kink pattern, cf. Fig. 3. Finally, the middle one is more 
complicated but one can easily guess that it is a six-kink pattern. 



Therefore, whatever the reason may be, it is certain that there exist various types of 
kinks that occur locally at different places in these rotational spectra. In other words, 
they are for sure not smooth (or analytic) functions of the spin. 

IV. CONCLUSION 

We have shown that the molecular rotational spectrum is not necessarily a smooth 
function of spin according to the measurements of R{I) and P{I) as well as the error 
estimate of the measurements. This is quite surprising because of our theoretical under- 
standing that the spectrum has to be an analytic function of spin /, which follows from 
the equation (|l]) obtained by assuming the Born-Oppenheimer approximation. Since it is 
rather difficult to believe that this approximation is violated, the first thought occurred 
to us was that there should be some error in the input file to our Fortran code because 
even a single error in typing data would create 'artificial' kinks. It may also be worth 
reconfirming all data at different laboratories with different experimental setups. 

In the previous section, we have shown relatively simple examples but there are actually 
many more complex cases. Fig. 5 shows the filtered one-point formula applied to the 0-0 
band belonging to the B^T,^ — -^^S^ system of the ^^Cu2 molecule P]. Nevertheless, 
it is likely that this result can also be decomposed and classified into various patterns 
as we have done in the previous section. We believe that this way of understanding the 
experimental data is suggestive and that intensive (phenomenological) studies along this 
line will give a useful clue. The aim of the present work is to initiate such a studyQ. One 
day, we may then be able to find the real physical implication lying behind this 'new' 
phenomenon, although whether or not it contains something new is an open question. 



^The Fortran source code used in the present work may be distributed upon request. Send an 
email to khara@physik.tu-muenchen.de or glalazis@physik.tu-muenchen.de. 
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FIGURE CAPTIONS 



Fig. 1 (a) Seventh order one-point formula and (b) tlie corresponding filtered formula 
applied to the YN data § 

Fig. 2 (a) Seventh order one-point formula and (b) the corresponding filtered formula 
applied to the YD data 

Fig. 3 Model spectra that have two, three and four successive kinks and the irregularities 
that occur in the one-point formula 

Fig. 4 Comparison between the V = H even spin data of the YD molecule (cf. Fig. 2b) 
and a three- kink 'theory' 

Fig. 5 Seventh order filtered one-point formula applied to the ^^Cu2 data 
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Fig. la 
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Fig. lb 



Even Spin 



Odd Spin 



0.006 




0.006 



0.004 



0.002 



-0.006 



4 8 12 16 20 24 28 32 36 40 



-0.002 



-0.004 



-0.006 




4 8 12 16 20 24 28 32 36 40 




4 8 12 16 20 24 28 32 36 40 



0.03 



0.02 



0.01 




-0.01 



-0.02 



4 8 12 16 20 24 28 32 36 40 



Spin I 



Spin I 



Fig. 2a 
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Fig. 2b 
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